Hole-pair symmetry and excitations in the strong-coupling extended t — J z model 
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We analytically calculate the ground state pairing symmetry and excitation spectra of two holes 
doped into the half-filled t — t' — t" — J z model in the strong-coupling limit (J z >> \t\, \t'\, \t"\). 
For the t' — t" — J z model, there are regions of d-wave, s-wave, and p-wave symmetry. We find 
that the t — J z model maps in lowest order onto the t' — t" — J z model on the boundary between d 
and p symmetry, with a flat lower branch of the pair excitation spectrum. In higher order d-wave 
symmetry is selected; however, we predict that the addition of the appropriate t' and/or t" should 
drive the hole-pair symmetry to p-wave. We perturbatively construct an extended quasi-pair for the 
t — J z model. We compare with analytic calculations for a 2x2 plaquette and numerical work, and 
discuss implications for the experimentally relevant parameter regime. 
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Although the issue has not been completely resolved, 
a variety of experiments has indicated that the pair sym- 
metry in the hole-doped cuprate superconductors is ei- 
ther pure d x 2_ y 2 or has a strong d x 2_ y 2 component 
Theoretical and numerical studies of the two-dimensional 
Hubbard, t — J, and related models believed relevant to 
the high-T c compounds have also suggested d x 2_ y 2 pair- 
ing ||-|6| , and Hubbard and t — J models on 2x2 plaque- 
ttes have recently provided an intuitive picture of how 
such a pair symmetry might arise However, there 

are few rigorous theoretical results in this general area. 

Recent experimental work has indicated that a pseu- 
dogap with the same symmetry as the superconducting 
gap can persist above T c in underdoped cuprate super- 
conductors This, along with the short high-T c 
coherence length , is generally consistent with a strong- 
coupling picture, where pairs can preform at T > T c [ p^[ . 
Numerical work has in addition suggested that the t — J 
and t— J z models have many similar properties ||,|l3].|l4| , 
and that the t — J z model may hence provide a suitable 
starting point for understanding t — J behavior p"5| . In 
this Letter, we consider two holes doped into the half- 
filled t — t' — t" — J z model in the strong coupling limit 
(J z » \t\,\t'\,\t"\). We calculate the symmetry of the 
hole pair in the ground state as well as the pair excitation 
spectrum. We consider first the t' — J z model, and show 
how singlet pairs can be constructed from our solutions. 
We next discuss the t' — t" — J z model, and then the t — J z 
and t — t' — t"— J z models. For the t — J z model, we per- 
turbatively construct an extended quasi-pair. As a step 
towards exploring the range of validity of our approach, 
we compare with results for a 2x2 plaquette and numeri- 
cal studies. Lastly, we discuss implications of our results 
for the physically relevant parameter regime, including 
the question of the sufficiency of the t — t' — t" — J model 
for capturing high-T c behavior. 

Specifically, we consider the Hamiltonian 



H — Hq + Hi + H2 



(1) 



where 



" , [ n x,y n x+l,y + n x,y n x,y+l )}, (2) 



H l = (-*) T, x ,y,„ {(4,3/,0-Cx+l.J/^ + H.C.) 

+ {ct,y,^x,v+l,a + H.C.)}, (3) 

H 2 = (-*') T,x,y,<r {(4,3/,<rC*+l,l/+l,<7 + H.C.) 



+ (ci, V ,aCx + l.y-l,a+H.C.)}, (4) 



and 



Hz = (-*") Ex, M {(4,,,/,+2,,,a + H.C.) 

+ {Cx,y,aCx,y+2,<j+ H.C.)}. 



(5) 



Here, x and y denote the coordinates of an LxL lat- 
tice with periodic boundary conditions and even L, and 



a = ±1 (|, I) refers to electron spin. 



L x,y, 



) , enforcing the condition of no double occupancy. 

&x,y = V 2 ( n x,yA ~ n x,y,i) an d n x,y = n x,y,1 + n>x,y,l- 

We do not explicitly consider here the spin-flip part of 
the magnetic interaction 



H± = 



Ji 



,y S - 



x+l,y 



X ,y X J 



H.c.}, (6) 



.r,i/ 



where S xy = 
t'-t" 



c x,yA°x,y,l an d S x 



Hq — J z Yl x ,y i(^x,y^x+l,y + ^x,y^x,y+l) 



c x, y ,i c *,yA- ( The 
full t — t' — t" — J model is recovered when J± — J z .) 

At half filling each site is occupied by exactly one 
electron, and the doubly degenerate ground state of 
Ho is then that of a Neel antiferromagnet. We choose 
|$ a > to denote the state with electron spins <j(x,y) = 
(— l) x+ y and | $(, > to denote the state with a(x,y) = 
(—l) x+y+1 . We define the operator a XtV = Cx,yMx,y) wl th 
a(x,y) = (~l) x+y , and the operator b x>y = Cx^^x.y) 
with a{x,y) = (— l) x +v+ 1 Although our calculations and 
results are independent of the ordering convention cho- 
sen, we will denote for specificity 

l $ a >= (al ii ...a t ljL )...(al i2 ...a t li2 )(4,i-a t i,i)|0 >, (7) 
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with an analogous definition for |$(, >. 

We now dope the half-filled state |$ a > with two 
holes and consider the strong-coupling limit (J 2 >> 
\t\, \t"\). In this limit, there will be an energy cost 
of order J z if the two holes are not nearest neighbors 
(n.n.). Hence, to zeroth order, the (highly degenerate) 
two-hole ground state is spanned by the set of all n.n. 
hole pairs. We denote the state with a horizontal n.n. 
hole pair at sites (x, y) and [x + 1, y) as 



(8) 



and the state with a vertical n.n. hole pair at sites (a;, y) 
and {x, y + 1) as 



>= 



L x,y+ 



(9) 



The \h x _ y >'s and \v x _ y >'s provide a complete, orthonor- 
mal basis for the two- hole ground state of H . 

It costs an energy of order J z if one of the n.n. holes 
hops to a n.n. site through the hybridization matrix el- 
ement t. However, there is no energy cost for hops cor- 
responding to t' or t" , as long as the two holes remain 
nearest neighbors after the hop. Thus, to lowest order 
in 1/Jj, it is only necessary to diagonalize the Hamilto- 
nian H2 + H3 in the subspace spanned by the \h x _ y >'s 
and \v x>v >'s; i.e., it is only necessary to consider the 
t 1 —t" — J z model. We note that in this limit the t'—t"—J z 
model becomes isomorphic to the strong-coupling limit of 
the antiferromagnetic van Hove model of |lq] . 

We consider first the t' — J z model, involving only the 
H2 (diagonal) hopping term. Defining 



1 



L ^ 

x,y 



(10) 



and 



1 ■r — ^ 2Tvik x x 27Tikyy 

\v kx , ky >= ^2^ e L e L Kv>. ( n ) 

X; y 

with k x ,k y = 0, 1...L— 1, we obtain the lowest order 
wave functions 



±sgn(i')e~^ \v kx , ky >} (12) 



with energies 



= ±4|t'| 



sm 



/7rk x \ /7rfc„ 
— - sin — i 



V L 



L 



(13) 



Since < sm(irk x /L), sm(irk y / L) < 1, the minus sign 
gives the branch of lower energy. The lowest energy state 
|V>o >, with energy — 4|i'|, occurs when k x = L/2 and 
k y = L/2 (i.e., (71", ?r)). Rewriting in terms of the a x ^s 
and neglecting overall phase factors, one obtains 



\4 a) >= 



LV2 



T, x .y {-t) x+v {ax+i, y ax, v 



- sgn(i') }\® a >. (14) 

When t' > (sgn(i') = 1), the sum over hole pair oper- 
ators in Eq. |lj changes sign upon a 90 degree rotation 
around a lattice point, giving the pair d-wave symmetry 
(specifically, d x 2_ y 2 |§||). When t' < 0, there are no 
such sign changes, giving s-wave symmetry (specifically, 
extended- s |^||). 

If one adds to Eq. |lj the appropriately-phased pair op- 
erator for two holes doped into the ground state |$h >, 
one obtains for t' > the usual (unnormalized) n.n. sin- 
glet d x 2_ y 2 pair operator 

-{c x ,y,]C x . y + l.l - Cx,y,lC x , »+l,t)}) (15) 

with t' < giving the analogous singlet extended-s oper- 
ator. With different relative phases, one can also obtain 
772 = triplet pairs; because quantum spin fluctuations 
are not included in the t — J z model, the two cases cannot 
be differentiated at this level. 



FIG. 1. Hole pair symmetry in the strong-coupling limit 
of the t' — t" — J z model as a function of t' and t" . Here, 
"D" denotes d x 2_ y 2, "P" denotes p x or p y , and "S" denotes 
extended-s. 

For the more general t' — t" — J z model, one obtains 
the (unnormalized) wave functions 

iV'fci- > = e-^(4±')s x s y \h K ,k y > 



(2t")(sl-sl)±T x ,y \v k , > (16) 



with energies 



s y) ^ T x,y 1 



(17) 



-2t")(l-4 
where s x — sm.(wk x / L) , s y — sm.(irk y /L), and 

T x , y = 2{{t''f{sl-slf+A{t'fslsiy. (18) 

As a function of t' and t" , we find that the ground state 
symmetry of the pair is as shown in Fig. 1. The p-wave 
pair operators can be either p x 



2 



Py 



l5> 

x,y 



+1,1/ 



-i,yJ 



*x 7 y— 1 j 



(19) 



(20) 



The p a states have energies independent of fc y , and the 
p v states have energies independent of k x . Both p-wave 
pair operators change sign under a 180 degree rotation. 

We next consider the strong-coupling limit of the t—J z 
model. To lowest order, we find that this maps onto the 
above strong-coupling limit of the t' — t" — J z model with 



a -t" -I t 

teff. - I, jr. - 3 I j 



(21) 



From Eq. [l?], the lower band of the pair excitation spec- 
trum then becomes flat, with wave functions 



> 



e L s y \h kx ,k y > 



e L s 



(22) 



Flat bands were also found 17 lq] for related models 
and/or treatments. In |l5) , a five- fold degeneracy of 
strong-coupling t— J z pairs of d or p symmetry was noted. 

We see from Eq. |2l] that, to lowest order, the strong- 
coupling t — J z model lies on the (rightmost) boundary 
in Fig. 1 between d-wave and p-wave symmetry. In the 
next higher order, neglecting constant additive terms, the 
energies of the lower band separate into 



& x ' & 'U 



&)(Jf)( 2 ~ C -- C 
{c 2 x + c y + 4c x Cy - 31c x 



31c y + 56}, (23) 



where here c x = cos(2irk x / L) and c y = cos(2%k y / L) . We 
then find (in agreement with p5|) that the pure d-wave 
(f > 0) state of Eq. |lj is selected as the ground state. 
However, the closeness to p-wave symmetry may provide 
an explanation for the low-energy p-wave "quasi-pair" 
peaks seen numerically in small t — J and t—J z clusters 
p3| . Because of this similar t — J and t — J z behavior, 
referring to Fig. 1 and assuming ground state pairs of 
pure symmetry, we predict that adding the appropriate 
t' and/or t" to the t — J z ort — J models with J z /t or J/t 
sufficiently large should drive the models to p-wave hole 
pair symmetry, and perhaps even p-wave superconduc- 
tivity. (In one dimension, a n.n.n. t' > will also give 
p-wave hole-pair symmetry in the J z » \t\,\t'\ limit.) 

One can perturbatively construct increasingly ex- 
tended quasi-pair states for the t — J z model. Combin- 
ing results for the n.n. d-wave pair operators for ground 
states |$ Q > and |$b >, one finds the lowest order cor- 
rection for the singlet pair operator of Eq. fa 



(-i)(x)(i)S> 



x,y,a 



{ ( C x+l 5 j/,(7 C ^+l ) y ) -< 7 C x,y+l,-a C x>y+l>v) 
Cx + 1 , y + 1 , a C-x , y , a 

(px+l,y,<7 C x-\-l,y,-<J ~ O x y _i _ a C x ^y-\^ 
Cx+l,y— l,crCx,y,(T } ■ 



(24) 



When operating on the appropriate Neel state, each of 
the above terms consists of a diagonal hole pair "dressed" 
with a singlet pair of electrons straddling the bond con- 
necting the pair of holes, as was recently found in numer- 
ical t — J simulations || . We note that the contribution 
from pairs a distance of two lattice sites apart, nominally 
also of order t/J z , vanishes identically in this order. This 
may provide an explanation for why only n.n. and diag- 
onal hole correlations appear to dominate in the t — J 
model near half filling for moderate to large J/t [|s|]l9f| . 

If one adds the necessary terms to the operator of Eq. 
pi] to impose rotational invariance, one obtains the com- 
posite pair operator invented in |lS) ] to give a diagonal 
singlet pair with d x 2_ y 2 symmetry. The non-invariant 
operator of Eq. which emerges naturally from per- 
turbation theory, also has d x 2_ y 2 symmetry. 

We also note that, since we calculate energy spectra 
and wave functions, our results and approach can be used 
to calculate finite-temperature and real frequency prop- 
erties. However, we do not pursue that here. 
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FIG. 2. Qualitative diagram of predicted hole pair symme- 
try for the t-t' - J z model. "D" , "P" , and "S" denote the 
same as in Fig. 1, and no prediction is made for region "IT". 

As one step towards investigating the range of validity 
of our approach, we performed analytic calculations of 
the t — t' — J z — J± model (see Eq. |)ona 2x2 plaque- 
tte. In general, we found that ground states remained 
smoothly connected as J z was reduced from strong cou- 
pling and also as J± was turned on. t — J numerical 
results |ll|[ljj also indicate that features of the strong- 
coupling limit may persist down to intermediate coupling 
(J/t sa 0.4—0.5). Together, the above support the strong- 
coupling t — t' — t" — J z model as a useful starting point 
for exploring the intermediate-coupling t — t' — t" — J 
model. Based on this and our strong-coupling results 
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(and, again, assuming pairs of pure symmetry), we show 
in Fig. 2 qualitative predictions of the hole pair symme- 
try for the t — t' — Jg model. We believe these predictions 
apply to the t — if — J model as well, with a compara- 
tively smaller p-wave region due to larger energy differ- 
ences between t — J p-wave and d-wave pair states [ [l3| . 
An additional t" > would enlarge the p-wave region. 

Reductions from CuC>2 three-band models , as well 
as comparison with ARPES results for a single doped 
hole (2J, suggest that \t'/t\ and \t"/t\ may be substan- 
tial, t > 0, and estimates for t' and t" are typically in 
the ranges if w (-0.1)* - (-0.5)*, and t" « 0^0 - (0.3)*. 
Both these signs of t' and t" could tend to drive the pair- 
ing symmetry to p-wave, raising the issue of the hole pair 
symmetry in the intermediate-coupling regime, (s-wave 
is also possible, though we believe it less likely at inter- 
mediate coupling.) It would be interesting to numerically 
explore whether the symmetry of two doped holes in the 
t — if — if 1 — J model is in fact d-wave for the experimen- 
tally relevant values of t,t',t", and J (e.g., J/t « 0.4). 
Drawing conclusions from exact diagonalization may be 
challenging due to finite-size effects: either s-wave or p- 
wave hole pair symmetries were found on lattices of sizes 
16, 18, and 20 for the t — if — J model with realistic 
parameters [ p2| . Another possible tool is higher order 
numerical ground state perturbation theory |lq| . If the 
symmetry were established to be p-wave rather than d- 
wave, it would suggest that the t — t' — t" — J model by 
itself could be incomplete as a model for high-T c super- 
conductivity. In that case, one possibility for restoring d- 
wave symmetry could be the addition of electron-phonon 
coupling in the d-channel |p3| . In either case, it may 
also be of interest to explore whether the existence of or 
nearness to p-wave symmetry, which effectively reduces 
the dimensionality of the hole pair wave function from 
2D to ID, might play a role in the "striping'' recently 
observed in certain of the high-T c cuprates |24j . 

In summary, we have investigated analytically the pair 
symmetry and excitation spectra of two holes doped into 
the half-filled t — t' — t" — J z model in the strong-coupling 
limit. In lowest order, this reduces to considering the 
if — t" — J z model, where we found regions of d-wave, 
s-wave, and p-wave symmetry We next found that the 
t— J z model in lowest order was on the boundary between 
d-wave and p-wave pair symmetry, with a flat lower pair 
excitation spectrum. In higher order, d-wave pairing was 
selected. However, because of the closeness to p-wave 
symmetry, we predict that the appropriate if and/or t" 
added to the t — J z or t — J models with intermediate to 
large J z or J should drive them into p-wave pairing, and 
perhaps even p-wave superconductivity. We constructed 
a perturbative correction to the nearest neighbor d-wave 
pair, and compared with the d-wave composite opera- 
tor invented in |19|. We explored ranges of validity of 
this perturbative approach using a 2x2 plaquette and 
results from other work Jl5| , ^9[ . Lastly, we discussed 
implications for the experimentally relevant parameter 
regime. These included the possibility of p-wave symme- 



try for two doped holes, which would suggest that the 
t — t' — t" — J model could be incomplete as a high-T c 
model, and the possible relevance of our results to "strip- 
ing" . 
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